Sign of the crossed conductances at a FSF double interface 
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Crossed conductance in hybrid Ferromagnet / Superconductor / Ferromagnet (FSF) structures 
results from the competition between normal transmission and Andreev reflection channels. Crossed 
Andreev reflection (CAR) and elastic cotunneling (EC) between the ferromagnets are dressed by 
local Andreev reflections, which play an important role for transparent enough interfaces and in- 
termediate spin polarizations. This modifies the simple result previously obtained at lowest order, 
and can explain the sign of the crossed resistances in a recent experiment [D. Beckmann et a/., 
cond-mat/0404360 . This holds both in the multiterminal hybrid structure model (where phase 
averaging over the Fermi oscillations is introduced "by hand" within the approximation of a single 
non local process) and for infinite planar interfaces (where phase averaging naturally results in the 
microscopic solution with multiple non local processes). 
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I. INTRODUCTION 



Andreev reflection^ is the mechanism by which charge 
is transported at normal metal / superconductor (NS) 
interfaces at voltages below the superconducting gap. 
A spin-up electron from the N electrode is reflected 
as a spin-down hole while a Cooper pair is transfered 
in the superconductor. Multiterminal geometries are 
of particular interest. For instance the sign of the 
crossed conductance between two electrodes connected 
to a same superconductor can be either positive or neg- 
ative depending on the relative strengths of the "nor- 
mal" (electron-electron) and " anomalous" (electron- hole) 
channels 3 ^* 5 -. Multiprobe structures are especially in- 
teresting when the electrodes are ferromagnets. Actu- 
ally, Andreev reflection also exists in ferromagnet / su- 
perconductor (FS) junctions but is suppressed as the 
spin polarization of the ferromagnet increases^* 7 ^. It 
was shown 9 ! 10 ! 11 . 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 19 ! 20 ! 21 that "crossed" An- 
dreev reflections (CAR) between two ferromagnetic elec- 
trodes, where spin-up electron and outgoing spin-down 
hole belong to different electrodes, give rise to interest- 
ing magnetoresistance effects. For instance, the latter 
process is favored by antiparallel spin polarizations. The 
fundamental problem of calculating the scattering matrix 
Sij of a NSN structure has been solved in some limits. 
The authors of Ref. |2(] use a "fork" geometry where 
the two electrodes merge at the superconductor. In con- 
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trast, in Rcfs. IlIollTll (like in a recent experiments!) , 
two distinct F/S contacts (denoted here as a, b) are sep- 
arated by a distance R. Yamashita et ali^S used a gen- 
eralization of the BTK approach in this geometrji^ 3 .. As 
a common result, CAR is possible if R < £, the super- 
conducting coherence length, but is also reduced by an 
algebraic factor. Andreev reflection becomes in this ge- 
ometry a genuine mesoscopic effect, non-local at the scale 
of the physical contacts. On the other hand, in the nor- 
mal channel, transmission between the two electrodes is 
also possible. Since it is spin-conserving, it is favored by 
parallel spin polarizations. In a tunnel model for the two 
F/S junctions, CAR and EC involve the virtual creation 
and propagation of a quasiparticle in the superconduc- 
tor. For the normal channel, this is similar to cotunnel- 
ing, introduced in the context of transmission across a 
Coulomb blockaded quantum doli^. The superconduct- 
ing gap replacing the Coulomb energy, and the normal 
crossed process has been called "elastic cotunneling" in 
Ref.E 

Other manifestations of spatially separated pair corre- 
lations were obtained in the study of equilibrium prop- 
erties of FSF trilayera 11 ! 25 ! 2 "! 2 ?! 28 ! 29 . It was shown 11 ! 25 ! 2 " 
within a model of multiterminal hybrid structure that 
the self-consistent superconducting gap can be larger in 
the parallel alignment. The same result was obtained 
for the FSF trilayer with atomic thickness, for half- metal 
ferromagnets^ 7 . and Stoner ferromagnets 2 ^. However sim- 
ulations with a finite thickness 2 ^ 9 , showed that pair- 
breaking dominates for strong ferromagnets as the thick- 
ness of the superconductor is larger than the Fermi wave- 
length (with therefore the superconducting gap larger in 
the antiparallel alignment). 

The goal of the theory of crossed conductances at a 
FSF double interface is to calculate the conductances of 
a FSF structure for any relative alignments of the spin 
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polarizations. In a first approximation we neglect the 
dependence of the superconducting gap on the relative 
spin orientation in the ferromagnets, which is a con- 
sistent assumption if the size of the contacts is much 
smaller than the superconducting coherence length^. 
The first approach to this problem was through Lan- 
dauer formalism^. Lowest order perturbation theory^ 
gives an interpretation in terms of CAR and EC pro- 
cesses. The effect of non collinear ferromagnets was also 
investigated", with the aim of describing transport 
of Cooper pairs at the interface between a superconduc- 
tor and a ferromagnet containing a domain wall. The 
Josephson effect between two superconductors connected 
by two spatially separated conduction channels was also 
examined^ and it was found that there is no Joseph- 
son effect within lowest order perturbation theory unless 
the length of the ferromagnets is smaller than the elastic 
mean free path, a condition that is not usually verified in 
experiments. Toy models for more complicated geome- 
tries involving Aharonov-Bohm effects related to crossed 
correlations were also investigated 15 . Disorder effects 
were also discussed recently for tunnel interfaces" and 
it was found that the geometrical reduction is less severe 
in the presence of disorder in the superconductor (dirty 
limit). Another geometry with a normal metal island 
connected to one su per conductor and two ferromagnets 
was proposed in Ref.[jjJ Noise correlations were also dis- 
cussed^ and recently within lowest order perturbation 
theory^. 



II. PRELIMINARIES 

In this section we provide the form of the Green's func- 
tions and Hamiltonians that we use throughout the arti- 
cle. 



A. Hamiltonians 

The superconductor is described by the BCS 
Hamiltonian^i: 

Hbcs = 51 e ( fc ) c k,^ Ck ' ff (!) 

+ A H ( c M<i + C k,lCk, T ) , 
k 

where e(k) = h 2 k 2 /2m is the free electron dispersion re- 
lation and A the superconducting gap. The ferromagnets 
are described by the Stoner model 

Wstoner = ^ ^tyc^C^ (2) 
k,(T 

- h ex ^2 (c£ T c k)T - c+ a c ka ) , 

k 

where h cx is the exchange field. 



CAR and EC were probed in a recent experiment by 
Beckmann et alm^. Driving a current through one con- 
tact induces a voltage in the other one. The overall ex- 
perimental results (the non-local resistance) are in agree- 
ment with theory except for the sign of the effect: it was 
predicted theoretically- 'SS, that for tunnel barriers the 
induced currents in the parallel and antiparallel spin ori- 
entations have an opposite sign whereas in experiments 21 
they have the same sign. Here we resolve this apparent 
contradiction by noting that the large interface trans- 
parencies used in experiments imply that the CAR and 
EC processes are "dressed" by local Andreev reflections 
at the two FS interfaces. This dressing is easily tran- 
scripted at zero temperature in terms of a perturbative 
expansion for the Keldysh Green's functions. Due to the 
strong damping of quasiparticles in the superconductor, 
a reasonable approximation involves a single non local 
propagator while local ones are treated to all orders. Wc 
also carry out numerical simulations of infinite planar in- 
terfaces in which multiple non local processes are taken 
into account and find a good agreement with the analyt- 
ical approach. 

The article is organized as follows. Preliminaries are 
given in section [H] The properties of multitcrminal hy- 
brid structures are investigated in section II I II Infinite 
planar interfaces are investigated in section IIVI Con- 
cluding remarks are given in section Ivl 



B. Green's functions in reciprocal space 

For collinear magnetizations the Nambu representa- 
tion corresponds to two 2x2 matrices, one in the sector 
S z = 1/2 and the other in the sector S z = —1/2 (S z is the 
projection of the spin along the quantization axis chosen 
parallel to the exchange field). The advanced Green's 
function in the sector S z = 1/2 is given at zero temper- 
ature by 

gA y ( t ,t') = -ie(t-t') (3) 

(K T (*),C+ T (t')}) ({Cx, T W, C y,l(*')}> 

({<j(*),< T (*')}) ({<;(*). Cy,lV)}) 

where x and y are two arbitrary sites and {.,.} is an 
anticommutator. Using the Hamiltonian one obtains 
the advanced Green's function in reciprocal spaced 2 .: 
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FIG. 1: (a) Schematic representation of the Andreev reflection 
geometry in a FSF structure, (b) Schematic representation of 
the tight-binding Hamiltonian cubic lattice and notation for 
the sites a, a, f) and b. (b) corresponds to a cut in the (y, z) 
plane. Voltages V a and VJ, are applied on the two ferromagnets 
and a voltage Vs is applied on the superconductor. 



where ps is the chemical potential in the superconductor, 
£ = h 2 k 2 /2m—ps is the kinetic energy with respect to the 
Fermi level, Ek = \J A 2 + £| is the quasiparticle energy, 



U% = (1 + £k/E)/2 and v 2 k = (1 - S, k /E)/2 are the BCS 
coherence factors, and 775 is a small parameter related to 
inelastic processes in the superconductor (a typical22iM 
value in experiments is r)s = 10~ 2 A). 

The "11" component of the Green's function of a fer- 
romagnet is given by 
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(w — pL a ) — £ + h c 



irjF 



(6) 



where fj, a is the chemical potential in the ferromagnet 
"a". A similar expression is obtained for g 2 ' 2 . The pa- 
rameter r)p contains information about decoherence in 
the ferromagneli^. 



C. Green's functions in real space 

It is useful in transport theory to manipulate Green's 
functions parameterized by spatial variables. To carry 
out the Fourier transform we define a short-scale cut-off 
corresponding to the lattice spacing ao of a cubic lattice 
tight-binding model where electrons hop between neigh- 
boring discrete sites like in Ref. |3^ (see Fig. [IJ. The 
tight-binding and continuous models are equivalent since 
we consider low energies compared to the band- width. 

The local Green's functions of a superconductor is de- 
fined as 



d 3 k 



(7) 



where a is a site in the superconductor (see Fig.^), and 
where g(£,u>) is given by Eqs. Q and (J5J. Evaluating 



the integral over wave vector leads to 



(8) 



y/A 2 - (w - fisf 



-(cj-hs) A 
A -(uj-fis) 



where ps is the normal state density of state, u> is the 
energy and A is the superconducting gap. The Green's 
functions of the superconductor take the formal 



^ (w) = fcw? exp 

sin (kpR) 



R 




(9) 



where a and are two sites in the superconductor (see 

Fig. HJ). CM = hv F /^/A 2 - (lu- ps) 2 is the BCS co- 
herence length. The local Green's functions of the ferro- 
magnet is given by 









Pa,l 



(10) 



where p a ,\ and p a ,i correspond to the spin-up and spin- 
down density of states in the ferromagnet "a" . A spin 
polarization P corresponds to p a ,\ — pf(1 + Pa) and 
Pa, 1 = Pf(1 — Pa)- We discard the energy dependence in 
(|10fl since we are interested in energies of order A, much 
smaller than the exchange field /i ox related to the spin 
polarization by P ~ h cx /ep, for h cx small compared to 
t F . 



D. Mixed Green's functions 

To describe infinite planar contacts we use Green's 
functions parameterized by the distance R along the z 
axis and the two wave- vectors k x and k y . This parame- 
terization is well suited for a situation where translation 
invariance holds in the [x, y) plane parallel to the in- 
terface but not in the z direction perpendicular to the 
interface. The kinetic energy £ is separated into the sum 
of f|| = h 2 (k 2 x + k 2 y )/2m - p s and £ x = h 2 k 2 /2m. The 
mixed Green's function is defined by 



9(R,t\\,u) = 



J ^e^g(^u), 



(11) 



where g(£,u) is defined by Eqs. and (jSJ. After evalu- 
ating the Fourier transform with respect to k z we obtain 



g A {R,t h u)= V^Mexp^fl) 



\ji\J A 2 - (u - us - ivs) 2 - T\\ 



(12) 
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with 



'-^i^(uj-[is-ivs) 2 ~A 2 -^. (13) 



by 



The mixed Green's function of a ferromagnet is given 



\j2ma\ 
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(14) 



^J-(uj - p a ) + ir)F +£\\-h' 



with 



K = 



'■^Juj ~Ha~ ir/F + h e 



(15) 



E. Transport properties 

The Green's functions Gij of the connected system are 
obtained by solving the Dyson equation. In a compact 
notation the Dyson equation takes the form G — g + 
g (g> £ £g) G, where £ is the self-energy corresponding to 
the tunnel Hamiltonian, ® is a summation over spatial 
variables and a convolution over time variables, and g 
is the Green's function of the disconnected system with 
£ = 0. 

Tranport properties'^ are obtained by evaluating the 
Keldysh Green's function 



[i + G 1 



(16) 



The spin-up current through the link a-a is given by 



Ia ' a ~ 2h 

+ {h e 



dwTr j t a , a G+; Q (uj) - f a , G+' Q (uj) <t z | 

■* -hex), (17) 



where the trace is a summation over the "11" and "22" 
components of the Green's function in the Nambu repre- 
sentation, and <7 Z is one of the Pauli matrices. The term 
(h C x — * -ft C x) corresponds to a summation over the "33" 
and "44" components in the 4x4 Nambu representation 13 
in the sector S z — —1/2. For a normal metal the four 
components of the current are equal and we recover the 
usual form of the transport formula'^. 



III. MULTITERMINAL HYBRID STRUCTURES 

A. Transport formula 

We start with the geometry on Fig.^ In the Hamilto- 
nian approach, the exact geometry is unimportant : the 



results are the same for two lateral contacts (like in the 
experiment^) or opposite contacts (Fig. P). Moreover 
we suppose in this section that the distance between the 
contacts is large compared to their transverse dimension 
so that the distance between sites a and (3 (see Fig. ^) 
can be considered as approximately independent on the 
choice of the sites a and /3 at the interfaces. This hy- 
pothesis does not hold in the infinite planar geometry 
(see section ITVT) . The total current flowing between the 
two electrodes is the sum of two contributions : i) the 
CAR current (Andreev process) circulating from both 
electrodes " a" and " b" , and ii) the EC current flowing 
from " a" to " b" . Those are given by 



ICAR.(V a ,Vb) = 2TT 2 tltl^-p a ^p bA 



(18) 




2,1,-R/ 



and 
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Eq. i|17|) is valid for all values of the voltages V a , V b and 
Vs. Eqs. (|18fl and (|19|l arc valid for all values of V a and V b 
and we supposed that Vs — since one reference voltage 
can be chosen equal to zero. 

The crossed conductance is defined^ as Q a b(V ai V b ) = 
dI a , a (V a ,V b )/dV b , where I a , a (V a ,V b ) = I C AR(V a ,V b ) + 
lEc(V a ,Vb) is the total crossed current through the link 
a-a. In the following we derive the currents (|18|l and Ijl9(l 
with respect to the voltage V b . One contribution to the 
derivative is due to the upper and lower bounds in the in- 
tegrals. A second contribution is due to the dependence 
of the Green's functions on V b . However we deduce from 
Eq. I© that g b ^ b is independent on V b because the chemi- 
cal potential is subtracted both in the energy uj — p,b and 
in the kinetic energy £. This second contribution to the 
derivative is thus vanishingly small. The crossed conduc- 
tance can be expressed in terms of the spin polarizations 
P a and Pb: 



Gab = ATT 2 t 2 a tl-p Fa PFb{ 

h 



(20) 
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FIG. 2: Variation of the normalized linear crossed conduc- 
tances XQp(P) (solid line) and XQap(P) (dashed line) as a 
function of the spin polarization P, for r = 0.25 (a), r = 0.5 
(b), r = 0.75 (c), r = 2 (d). Note the different scaling factors 
on the conductance axis, and r = 1 



+(eV b -> -eV b )} . 

Averaging over the Fermi oscillationsiSii^ is used to sim- 
ulate an extended contact as the sum of channels in par- 
allel with a distribution of length R: 



r,l,2,A, \^,2,l.i?./ \ 



1 

~ar 



R+AR/2 



dR 



R-AR/2 



(21) 
(R,cj), 



with AR = X F = 2n/k F . 

The tunneling limit involves bare Green's functions, 

, . , ,. r 1.2. A, \ 2.1,R, \ 1,1. A i \ 1.1, R, ~\ 

which satisfy g a0 (uj)g p \ a M = g^J {uj)g p \ a (w), 
therefore the crossed conductances for parallel (P a = Pb) 
and antiparallcl (P a = —Pb) spin polarizations are op- 
posite. As a result, the crossed conductance is zero in 
absence of spin polarization^. Similarly, the crossed con- 
ductance at a mixed FSN interface (where only one elec- 
trode is ferromagnetic) is zero. As shown below, these 
symmetry properties do not hold anymore for transpar- 
ent interfaces, where dressing of the quasiparticle propa- 
gators occurs in the superconductor. 



B. Perturbative expansion in l/(kpR) 

Now we revisit the perturbative expansioni2*i2iiii used 
to describe crossed Andreev reflection. We suppose that 
the contacts a-a and 6-/3 can be highly transparent so 
that multiple Andreev reflections take place locally. How- 
ever hp R is large so that in a first approximation a single 
non local process is included in G a ^. This leads to 
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FIG. 3: (a) A CAR process to order r4. (b) A CAR process to 
order t 8 involving two local Andreev reflections. The Nambu 
indices are indicated at each vertex. A similar diagram can 
be drawn for EC process. 



where Mai (iV?]) describes the dressing by multiple 



Andreev reflections at the contact a-a (6-/3). Ma) a is de- 
termined by the equation Ga}a = Ma) a g a ^ a , where Ga} a 
is the fully dressed propagator at site a with % $ = 0. 
We find easily 



K°) - 



N, 



0.0 



9oi-,ct^a,a9a,a^a,a. 



(23) 
(24) 



Generalizing Ref. [jjj we evaluate the phase averaging of 
the following Green's functions: 



1,1 2,2 
9a,09 a ,0 

171 O" 
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(25) 
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■ exp 



■ exp 
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— (26) 
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(27) 



£{u))J A 2 



-loA 



We use Eqs. P5 )) -l|37 )l and Eq. (J22J| to evaluate the phase 
averaged Green's function in the CAR and EC currents 
given by (JTSJ and JT^J. 



C. Crossed conductance in a FSF structure 

1. Normalization of the tunnel amplitudes 

For a NN contact, highly transparent interfaces corre- 
spond to t — to, where to is such that ir 2 t 2 p 2 N — 1. The 
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NN contact conductance is given by Gnn = (2e 2 /h)a, 
with2i 



x=0.25 



t=0.5 



a 



4(*/*o) 



(i + (*Ao) 2 ) 2 



(28) 



The conductance of a NF contact is given by Gnf 
(e 2 //i)(a| +cxj.), with 



aj = 
"1 = 



47r 2 t 2 PA rp F (l + P) 

[l + ^ 2 t 2 pwPF(l + -P)] 2 

47r 2 t 2 Pjv p F (l-P) 
[l + 7r2t2 pAr p F (l-P)]2 



(29) 
(30) 



The value to of t corresponding to a perfect transmis- 
sion in the spin-up channel is given by 



ir 2 t 2 (P)p NPF {l + P) = 1 



(31) 



The conductance of a SF contact with perfect transmis- 
sion in the spin- up channel (t = t (P)) is given by 



Gnf = t" 1 



1 



1=p\ 
i+P J 



(32) 



varying between e 2 /h for half-metal ferromagnets and 
2e 2 /h in the absence of spin polarization. In the following 
we normalize the hopping amplitude t to the maximal 
value to(l) of t (P): t = ri (l), with r between and 1. 



2. Linear crossed conductance 

Taking u> <C A, the linear crossed conductances corre- 
sponding to Vb = are given by 



A<?ap 
XGp 



r 4 (l-P 2 ) 
V 2 (P) 



e 2 t 4 
/i £> 2 (P) 
,e 2 r 4 



/i £> 2 (P) 
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P 2 



(33) 
(34) 



P 2 (P) 
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with £>(P) = 1 + t 4 (1 - P 2 )/4, A = 2(k F R) 2 exp (2P/£ ) 
in the ballistic limit, and A = 2(kpl e )(kpR) exp (2i?/£o) 
in the diffusive limi^, where l e is the elastic mean free 
path and £o the coherence length at zero energy. The 
variations of XGap (P) and XGp (P) are shown on Fig. El 
for increasing values of t. First, one sees that at P = 
the crossed conductance is negative, as if a normal metal 
would replace the superconductor. By contrast, in the 
tunneling approach, an exact symmetry holds between 
CAR and EC, so that the crossed conductance is zero for 
P a = or Pi, = 10 ' 12 . The observed trend is more ap- 
parent for transparent contacts, and means that EC pro- 
cesses dominate over CAR. We obtain Gap(P) — —Gp(P) 
only for small values of r (r = 0.25 on Fig. EJ-(a)). As r 
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FIG. 4: Variation of the normalized crossed conductance 
X(Vb)G(Vb) = XdI a /dVb(Va = 0,14) as a function of the nor- 
malized voltage Vb/ A, for P = 0.5 and increasing values of the 
dimensionless tunnel amplitude r. The parallel (antiparallel) 
alignment corresponds to the solid (dashed) line. 



is increased Gap{P) changes sign as shown on Figs. 
(b), 13(c) and Gl-(d). For half-metal ferromagnets we 
have Gap(P) = —Gp(P) for all values of r, as it can 
be seen from Eqs. I)33|) and (|34|) . The key role in the sign 
changes is played by processes of order r 8 like the one 
represented on Fig. 13(b). Due to the opposite signs of 

(da^) 2 and g^pQa^ forw = ( see E q s - local 

X 2 j4. 2 1 

Andreev dressing effects tend to decrease G a '^ Gg'^ 

and increase G^'^G^'^, therefore CAR is weakened 
and EC reinforced. 



3. Crossed conductance versus voltage 

The variations of the crossed conductances as a func- 
tion of the voltage Vb applied on electrode "b" are shown 
on Fig. 2] for P = 0.5. For small interface transparen- 
cies (r = 0.25 and r = 0.5) the crossed conductances 
in the parallel and antiparallel alignments are approxi- 
mately opposite in the entire voltage range. The crossed 
conductance at Vb = A is vanishingly small, both in the 
parallel and antiparallel alignments and for arbitrary spin 
polarizations. 

For half-metal ferromagnets the crossed conductances 
take the form 



X(V b )GAp(V b ) 



-X(V b )Gp(V b ) 
A 2 



v b 2 



v b 2 



A 2 -If 



(35) 

-2 



where A (Vb) has the same expression as A except that 
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£o is replac ed by the coherence length at a finite en- 
ergy £oA/a/ A 2 — V b 2 . Eq. (|3*5f is approximately equal 
to (8/r 4 )(l - Vb/A) if A - v < r 4 A, and approximately 
equal to 4r 4 A 2 /(A 2 - V b 2 ) if A - V b > r 4 A. For r small 
there is thus a maximum in the crossed conductance^ at 
a voltage V& ~ A(l — t 4 /2). This argument for P = 1 is 
compatible with the behavior of the crossed conductances 
for V b < A and P = 0.5 on Fig. H 



^. Crossed conductance of a FSN or NSF junction 

The linear crossed conductance of a FSN (NSF) junc- 
tion 



t/A=3.2, k F R=10 



t/e F =0.13, k F R=10 J 



Sr 8 (l - P 2 



A (l + r 4 (l-P 2 )/4)(l + r 4 /4) 



(36) 



is proportional to r 8 and thus very small for tunnel inter- 
faces. However it can take measurable values for highly 
transparent interfaces and with a weakly polarized ferro- 
magnet or even a normal metal. 



D. Interpretation of the dressing by local processes 



One of the terms contributing to G a ' ^ a 
parallel alignment is 



the 



M (o),i,2 2,1 ^(oj.i.ijif (0),i,2 2,1 ^(0,1,1 f371 
from what we deduce one of the terms of order i 4 t 4 con- 



tributing t 2 a t 2 bPaAPb , T G£ Jj A Gj£ A ■ 



2,2 2, 2, A .2,2 2,1 
:,ai/a,a a.a9 a.J3 



t 1 ' 1 o 1 ' 1 ^ 1 a xa t 

,1,1 1,1,1,1 1.2,2.2 2. 2, i? ,2.2 2,1 

Xt /3,bPb,b l b,l39f3,f3 T f],b9b,b t b,p9p t a 



(38) 
(39) 



This process corresponds to a diagram in which a spin-up 
electron from electrode "a" is locally Andreev reflected 
as a spin-down hole that makes an excursion in electrode 
"a" , is converted in a spin-up electron in electrode "b" 
through a CAR. The spin-up electron in electrode "b" 
undergoes a local Andreev reflection, is converted in a 
spin-down hole that is transformed in a spin-up electron 
by a CAR. This process is not possible for half-metal fer- 
romagnets, in which case the linear conductance is equal 
to the bare conductance, proportional to r 4 . 



IV. INFINITE PLANAR INTERFACES WITH A 
BULK SUPERCONDUCTOR 

Now we consider infinite planar 2D interfaces connect- 
ing a bulk 3D superconductor to bulk 3D ferromagnets. 
Within this approach we do not have to impose "by 
hand" phase averaging as we did in the previous section 
(see Eqs. J250 - H270 ). We use the mixed Green's functions 
introduced in section lTl Dl The evaluation of the integrals 
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FIG. 5: Variation of the total current I a through electrode 
"a" evaluated with eV a = and eVb = A/2 as a function of 
the exchange field, with A/e.F = 4 x 10 -3 . The parallel (an- 
tiparallel) alignment corresponds to the solid (dashed) line, 
(a), (b) and (c) correspond to a full Green's function calcula- 
tion with CAR and EC to all orders, (d) corresponds to the 
approximation given by Eq. I]22l . We use rjs/A — 10 _1 and 
t)f/£f = 0.4 but obtained similar results for ?;f/A = 10 -1 . 



in the expression of the CAR and EC currents is briefly 
described in Appendix^ A realistic model would involve 
an estimate of the inverse proximity effect, i.e. how the 
superconducting gap depends on the relative spin orien- 
tation of the ferromagnetic electrodes 2 ^S*SLSSi 2 2. Here 
we consider the infinite planar limit just as a test of the 
approach in the preceding section and do not impose self- 
consistency on the superconducting gap. 

We calculate the crossed current through electrode "a" 
with a finite voltage Vb on electrode "b" , and with V a = 
and suppose that the distance between the ferromagnets 
is much larger than the Fermi wave-length. The varia- 
tions of the current evaluated at Vb = A/2 as a function 
of h cx /A is shown on Fig.JSJfor different values of the dis- 
tance R between the contacts. For small values of tj A we 
obtain one Andreev bound state for weak ferromagnets 
for which h ex and A have the same order of magnitude 
(see Fig. |5}(a)). Andreev bound states for weak ferro- 
magnets were already discussed elsewhere 2 ^. For larger 
values of t/A and larger values of R/ao (see Fig.El-(b)) we 
obtain Jap — Ip up to relatively large values of h ex / A 
and Jap — —Ip for the largest values of h ex /A. The 
overall variation on Fig. EJ-(a), Fig. 0(b) and Fig. EJ(c) 
looks like Fig. EKa), Fig.H(b) and Fig.|3(c). Fig. 13(d) 
corresponds to the same parameters as Fig.[5J-(a) except 
that we used the approximation l|22|l of a single non lo- 
cal process on Fig.OD-(d) whereas Fig.GJ(a) was obtained 
without approximation. The two variations are almost 
identical, except for the regime h cx <C A. The agree- 
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ment is even better for larger values of R (not shown on 
Fig. [SJ. This illustrates the validity of the approxima- 
tion 1)22(1 in which we keep a single CAR. 



V. CONCLUSIONS 

To conclude we have proposed a model that may con- 
stitute an explanation to the signs of the crossed currents 
in the recent experiment by Beckmann et al. 21 . In the ex- 
periment, a current is injected through electrode 6, and a 
voltage V a is measured. Assuming that the local conduc- 
tances Gaa, Gbb are much larger than Gab, Gba, one simply 
gets the non-local resistance ~ Sab V /GaaGbb- 

For small interface transparencies we recover the results 
of perturbation theory to order t\t\ : the crossed con- 
ductance in the antiparallel alignment is opposite 
to the crossed conductance G\\ in the parallel alignment, 
and is zero if P a = or = 0. For larger interface trans- 
parencies, the propagators of crossed Andreev reflection 
and elastic cotunneling are dressed by local Andreev re- 
flections. As a result Gab can become negative as the 
interface transparencies are increased, whereas Q* b is al- 
ways negative. We also carried out simulations of infinite 
planar interfaces and found the same qualitative behavior 
without imposing phase averaging by hand and without 
using the approximation in which a single non local pro- 
cess is retained. 
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APPENDIX A: EVALUATION OF INTEGRALS 
FOR INFINITE PLANAR INTERFACES 

In this Appendix we give necessary technical details 
regarding changes of variable for evaluating the current 



with infinite planar interfaces. The "11" component of 
the CAR current takes the form 



r car .3,2 ,2 e / 2ma 2 \ 



eV a rD 

duj I df||(Al) 



-eV h 



-D 



and the "11" component of the EC current takes the form 



Pi.'iV - eK,£||)Pi',i(w _ eV b^\\) GiT' p (w,£||) 



Similar expressions are obtained for the "22" currents 
and for the currents in the spin-down sector. 

The density of state p°'°(w — eV^,^) and Pii(u> — 
eVb, deduced from Eq. ^] contain a square root sin- 
gularity. For instance 



n sju- eV a - §|| + h a 
After a change of variable takes the form 



du 



n 2 



dx 



(A4) 



^2x + u 2 



where explicit expressions of x(u>), x m in, Smax and Umax 
can be obtained in each case. For instance if — eVb < eV a 
we have x(ui) = —uj + e(V a + Vb)/2 — (h a + hb)/2 and 
u = ^Jlo — eV a — fii + h a . The integral l|A4() is then eval- 
uated by making the changes of variable v 2 = 2x and 
(u, v) = p(cos6>,sin#), therefore absorbing the square 
root singularity in a change of variable. 
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